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Abstract 

We use density functional theory to quantify finite size and shape effects 
for gold nanoclusters. We concentrate on the computation of binding energy 
as a function of bond length for icosahedral and cuboctohedral clusters. We 
find that the cuboctoheral gold clusters have lower energy for 13 atoms. For 
55 atoms we find that the icosahedral gold clusters have lower binding energy. 
We also introduce a one parameter family of geometries that interpolate be- 
tween the icosahedral and cuboctohedral clusters that is parametrized by an 
angle variable. We determine the binding energy dependence on shape as a 
function of the angle variable for 13 and 55 atom clusters with a minimum at 
the cuboctohedral point and icosahedral point respectively. We also compute 
the binding energy for the 147 atom gold nanocluster and show that the bind- 
ing energy of the icosahedral cluster is lower than the 147 atom cuboctohedral 
gold cluster. We also compute the binding energy of the A1155O2 molecule with 
possible applications to catalysis. 

1 Introduction 

An important application area for nanoscience are gold nanoclusters. These clusters 
can be used in energy science as catalysts and in medicine for drug delivery. It is 
also important to understand the basic physics and chemistry of gold nanoclusters 
to understand the limitations of the current applications as well as to suggest new 
ones. Experiments on gold nanoclusters have recently been performed on clusters 
with 102 gold atoms [1]. In addition increasing computational power has lead to the 
ability to simulate a nanocluster of hundreds of atoms on supercomputers running 
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scalable density functional codes. Thus this is an exciting time in the development 
of nanoscience where computation, theory, and experiment are helping to improve 
our understanding of these systems. 

The main tool that we use to computationally probe gold nanoclusters is density 
functional theory run on large parallel computers. Density functional theory is a 
computationally efficient approach to many body systems that maps the interacting 
many body Schrodinger equation to a single body noninteracting problem given by: 

h 2 \ 
- 2^ V 2 +V S {-?) J <M r>) = eiVi(^) 

where V s is defined by: 

V,(?) = V ext (T) + J yz^T^ + Vxc[n(T)} 

with Vexti'r') the external potential, Vxc[ n C?)] the exchange correlation potential 
and n(T^) the electron density. 

One typically uses a local density approximation and generalized gradient ap- 
proximation where the exchange functional is a local function of the density and 
gradient of the density. One then iteratively solves the one-body equation using 
the relation between the electron wave function and density to evaluate the the ex- 
change correlation potential. More details about the density functional approach 
can be found in many textbooks, for example [2]. 

Density functional theory has been implemented in several quantum chemistry 
codes. In this paper we use the NWChem quantum chemistry code distributed by 
PNNL [3]. We also use a large Blue Gene L parallel supercomputer that efficiently 
runs NWChem up to 2048 nodes or 4096 processors. A general description of the 
application of NWChem to materials science can be found in [4]. Outstanding fea- 
tures of the NWChem input file specify the geometry, basis, effective core potential, 
exchange and correlation functionals, DIRECT, SMEAR and MULT. We used the 
basis set specified by SBKJC VDZ ECP [5] and the pbe exhange and correlation 
functionals [6]. The DIRECT keyword indicates that all integrals are computed 
on the fly, the SMEAR keyword allows fractional occupation of molecular orbitals, 
and the MULT keyword defines the spin multiplicity [7]. In this paper we used the 
default smear value of .001 and multiplicity 2. 

Besides the practical application of gold nanoclusters little is known of the struc- 
ture of the energy landscape of theses systems even though the chemistry and mate- 
rial properties strongly depend on the position in this landscape. The main difficulty 
is that even with 13 atoms one has a 13 x 3 dimensional space to explore and this 
is difficult computationally as well as difficult to visualize. Thus the main purpose 
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of this paper is to find a simpler description of this landscape in terms of the size 
and shape of the gold nanoparticles for well known structures like the icosahedron 
and cuboctohedron. Previous density functional studies of gold nanoparticles are 
given by [8], [9], [10], [11], [12]. In particular [8] and [9] studied the relative energy 
between icosahedral and cuboctahedral gold clusters and our results are consistent 
with their results. The nanogold phase map in temperature and critical size has re- 
cently been determined based on relativistic ab initio calculation and verified using 
high resolution electron microscopy at finite temperature [13]. Their work indicates 
that different shapes of nanoparticle are more stable depending on the place in the 
phase diagram. By varying the size and shape of the gold naocluster we are able 
to determine the stable configurations as a function of critical size using density 
functional methods at zero temperature. 

For gold, relativistic effects are very important. In [14] density functional stud- 
ies indicate a qualitative difference between 55 atom gold and silver nanoclusters, 
and found gold nanoclusters have several low-lying low-symmetry structures and 
washed-out electron shell structure [14]. In [15] photo-electron spectroscopy and 
first principle calculations indicate a nonicosahedral disordered cluster with strong 
surface contractions due to relativistic effects. In this paper we use an exchange 
potential that includes relativistic effects. However we focus on two well known 
structures cubocahedral and icosahedral. This leaves aside the question of the true 
lowest energy state of the system, but allows us to study the relative stability of 
the two structures as well as the transition and intermediate geometries between the 
icosahedron and cuboctahedron. These type of Mackay transitions between different 
shapes have been studied for lead [16] and iron [17] and we extend this analysis to 
gold nanoclusters here. 

This paper is organized as follows. In section 2 we present the computation of 
the energy of 13 atom gold clusters in the shape of a cuboctahedron and icosahe- 
dron as well as intermediate structures. We derive a potentials which describes the 
dependence of the energy on the size and shape of the system. In section 3 we do 
the same thing for 55 atom gold clusters and compare with the previous section. We 
find that the cuboctahedron is more stable than the icosahedron for 13 atoms while 
the icosahedron is more stable than the cuboctahedron for 55 atoms. In section 4 
we study the binding energy of the AU55O2 molecule and applications to oxidation 
and catalysis. We calculate the potential associated with the separation of the two 
oxygens and show that the presence of the gold nanoparticle lowers the potential 
relative to the free 0%- In section 5 we state the main conclusions of the paper. 
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Table 1: Number of atoms N for a cuboctahedron or icosahedron for integer K 
which is the shell index. 



K 


N 


2 


13 


3 


55 


4 


147 


5 


309 




Figure 1: 13 gold atom nanocluster in the shape of an icosahedron (left) and cuboc- 
tahedron (right). 

2 Finite size and shape effects in 13 atom gold nan- 
oclusters 

In density functional theory (dft) one usually computes the energy E of the multiple 
atom system and derives the binding energy V from: 

V = E — NEatom 

where N is the the number of atoms and E a t om is the dft computation of the energy 
of a single gold atom. 

The total number of atoms in a cuboctohedron or icosahedron N is given by 
[18]: 

N = —K 3 — 5K 2 + —K — 1 
3 3 

It Table 1 we list the first few values of N for integer K which is the shell index. 

We show the 13 atom icosahedral and cuboctahedral geometry in Figure 1. For 
the case of K = 2 and N = 13 we list the coordinates for the cuboctahedron in table 
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Table 2: Coordinates for 13 gold nanoclusters with bond length a and shape of 
cuboctahedron, icosahedron and intermediate structure with shape parameter 9. 



cuboctohedron 



«73(0,1,1) 
a^(0,-l,l) 
a^(0,l,-l) 
a^(0,-l,-l) 

a^(l,l,0) 
a^(-l,l,0) 
073(1,-1,0) 

a 72(" 1 '- 1 '°) 
^(1,0,1) 

a^(l,0,-l) 

a^(-l,0,l) 

a^(0,0,0) 



icosahedron 



<4(0, 1, ^) 
a±(0, 

a±(0, 1, -99) 
a±(0,-l,-<p) 

4(i,^,o) 

ai(-l,v3,0) 
0^(1,-^,0) 
-<P,0) 
a±(ip,0, 1) 
05(^,0,-1) 
a±(-ip,0, 1) 
a±(-v?,0, -1) 
a|(0,0,0) 



intermediate 



a(0, cos 9, sin 0) 
a(0, — cos 9, s'm9) 
a(0, cos 0, — sin 6*) 
a(0, — cos 9, — sin9) 
a(cos 9 , sin , 0) 
a(— cos 0, sin 9, 0) 
a(cos 0, — sin9, 0) 
a(— cos 0, — sin 0, 0) 
a(sin#, 0, cos 9) 
a(sin 9, 0, — cos 6*) 
a(— sin 0,0, cos 6*) 
a(— sin 6*, 0, — cos#) 
a(0,0,0) 



Table 3: Data near the minimum of the binding potential for 13 atom gold nan- 
oclusters. 



shape 


bond length [nm] 


Binding Energy [eV] 


cuboctahedron 


0.276897 


-25.1064 


cuboctohedron 


0.280639 


-25.3946 


cuboctohedron 


0.28438 


-25.5349 


cuboctohedron 


0.288122 


-25.4893 


icosahedron 


0.285756 


-24.1683 


icosahedron 


0.291047 


-24.5989 


icosahedron 


0.296339 


-24.6706 


icosahedron 


0.301631 


-24.4062 



2 with bond length a. We list the coordinates for the 13 atom icosahedron also in ta- 
ble 2. Here ip = is the golden ratio which is approximately 1.618033988749895. 
For the icosahedron there are two bond lengths, a and approximately 0.951057a and 
we use the former. One can also consider intermediate shapes parametrized by a 
shape parameter 9 whose coordinates are given in Table 2. Here 9 = 45° corre- 
sponds to the cuboctohedron and 9 = arctani/? « 58.282525589° corresponds to the 
icosahedron. The list of coordinates for intermediate geometries is given in the third 
column of Table 2. 

In this paper we primarily study K = 2 and K = 3 which corresponds to 13 
and 55 gold atoms respectively. Higher values of K can also be studied and require 
large computational resources and scalable density functional codes. In section 3 we 
present some results for K = 4 and 147 gold atoms. 

Figure 2 shows the binding potential for the cuboctohedron and icosahedron 
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13 atom icosahedron 
13 atom cuboctahedron 



--X-- 



-X- 



0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 
Bond Length [nm] 



Figure 2: Binding energy as a function of bond length for 13 atom gold nanoclusters 
shaped like an icosahedron and cuboctahedron. 
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Figure 3: Binding energy as a function of the shape parameter 6 for a 13 atom gold 
cluster. The minimum occurs at 45 degrees which corresponds to the shape of a 
cuboctahedron. 
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Table 4: Coordinates for 55 gold nanoclusters with bond length a and shape of 
cuboctahedron, icosahedron and intermediate structure with shape parameter 9. 
Here n is the number of atoms in each class of coordinates. 



cuboctohedron 


icosahedron 


intermediate 


n 


13 atom cuboctohedron 
same as above scaled by 2 
a^(±2,0,0) 
a^(0,±2,0) 
a^(0,0,±2) 
a^-(±2,±l,±l) 
a^(±l,±2,±l) 
aj-(±l,±l,±2) 


13 atom icosahedron 
same as above scaled by 2 
ai(±29?,0,0) 
a±(0,±2<^,0) 
a±(0,0,±2v?) 
a±(±(l + <p),±tp,l) 
a±(±l,±(l + <p),±<p) 
a%(±<p,±l,±(l + <p)) 


13 atom intermediate 
same as above scaled by 2 
a(±2 sin 0,0,0) 
a(0,±2sin<9,0) 
a(0,0,±2sin6>) 
a(±(cos 6 + sin 9), ± sin 9, ± cos 9) 
a(± cos 6*, ±(cos# + sin 6*), ± sm.9) 
a(± sin 6*, ± cos 6*, ±(cos 9 + sin 6 1 )) 


13 
12 
2 
2 
2 
8 
8 
8 



structure. The difference in the potentials is rather small so we list the values near 
the bottom of the potential in Table 3. From Table 3 we see that the binding energy 
of the cuboctohedron shape for 13 gold atoms, -25.5349 eV, is greater in magnitude 
than that of the icosahedron, which was -24.6706 eV. Thus our calculations indi- 
cates that the cuboctahedron is more stable than the icosahedron for 13 gold atoms. 
The binding energy of the cuboctahedron shape minus the binding energy for the 
icosahedron shape for 13 gold atoms is AE = —.8643 eV. This is in agreement with 
[8] [9]. Figure 3 shows the dependence of the energy on the shape parameter 9. The 
potential shows a minimum at 9 = 45 degrees which corresponds to the cubocto- 
hedron. Thus the cuboctahedron is stable with respect to shape transformations 
relative to the icosahedron. Indeed from the figure the icosahedron is not a local 
minimum and is unstabe with respect to shape perturbations for the 13 atom gold 
nanocluster. 

3 Finite size and shape effects for 55 atom gold nan- 
oclusters 

One can proceed in a similar manner for the 55 atom gold cluster. The cuboc- 
tahedron and icosahedron 55 atom gold nanocluster are shown in Figure 4 and 
coordinates are given in Table 4. The binding energy a function of bond length is 
shown in Figure 5 and Table 5 for the cuboctahedron and icosahedron respectively. 
The data indicates that the icosahedron is lower in energy than the cuboctahedron 

with AE = E cuboctahedron ^icosahedron = 1.28019 eV. 

In order to further examine the stability of 55 atom gold clusters with respect to 
change in shape we plot the binding energy as a function of the shape parameter 9 in 
figure 8. The angle 9 is defined similar to the 13 atom case as being 45 degrees for the 
cuboctohedron, arctan ip for the icosahedron and intermediate values in between. A 
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Figure 4: 55 gold atom nanocluster in the shape of an icosahedron (left) and cuboc- 
tahedron (right). 
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Figure 5: Binding energy as a function of the bond length for a 55 atom gold cluster. 
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Figure 6: Binding energy as a function of the shape parameter 8 for a 55 atom gold 
cluster. The minimum occurs at 58.282525589 degrees which corresponds to the 
shape of a icosahedron. 



Table 5: Data near the minimum of the binding potential for 55 atom gold nan- 
oclusters. 



shape 


bond length [nm] 


Binding Energy [eV] 


cuboctahedron 


0.276897 


-129.961 


cuboctohedron 


0.288122 


-135.901 


cuboctohedron 


0.291864 


-135.319 


cuboctohedron 


0.299348 


-132.635 


icosahedron 


0.285756 


-132.143 


icosahedron 


0.291047 


-135.901 


icosahedron 


0.296339 


-137.181 


icosahedron 


0.301631 


-136.643 


icosahedron 


0.306923 


-134.659 



Table 6: Data near the minimum of the binding potential for 147 atom gold nan- 
oclusters. 



shape 


bond length [nm] 


Binding Energy [eV] 


cuboctahedron 


0.28438 


-383.226 


cuboctohedron 


0.288122 


-387.14 


cuboctohedron 


0.291864 


-387.988 


cuboctohedron 


0.295606 


-386.206 


icosahedron 


0.291047 


-385.21 


icosahedron 


0.296339 


-391.058 


icosahedron 


0.301631 


-391.033 
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simple prescription to define the shape parameter is to replace the occurrence of the 
golden ratio ip in the list of icosahedral coordinates with tan# and scale all the 55 
coordinates proportional to cos# to construct the interpolating nanostructures. The 
same prescription works also for the 13 atom C MS C (IS indicated in tables 2-4. The 
inner bond length for the 55 atom case is fixed at 0.951057ao where ao = 0.296339 
nm. The resulting plot indicates the the icosahedral structure is at lower energy 
than the cuboctahedral structure for nanoclusters consisting of 55 gold atoms with 
respect to change of shape. 

It should be emphasized however that we have just investigated the relative 
stability between two well known structures, the cuboctohedron and the icosahedron. 
There are many more structures that do not fit into these categories or even the 
interpolation between them. Some of these have lower energy [19]. Nevertheless the 
ability to computationally determine the potentials for size and shape of nanoclusters 
gives a first simplified look at the so called energy landscape of theses models. The 
demonstration of the ability to use efficient density functional codes and parallel 
computing to map out this landscape is one of the main results of this paper. 

One can extend the results of this paper to larger gold clusters. Table 6 shows 
the results of a binding energy calculations of 147 atom gold clusters. The table 
indicates that for the 147 gold atom case the icosahedron is lower in energy than 
the cuboctahedron with AE = E cuboctahedron - E icosahedron = 3.07 eV. 

4 Calculation of Au^Oi binding energy 

Besides pure gold clusters another area of interest is the use of gold clusters as 
oxygen catalysts involving the disassociation of O2 molecules. 

In [24] X-ray spectroscopy was used to show that 55 atom gold clusters have 
a maximum oxidation resistance and suggests that they may be effective oxidation 
catalysts. In [25] hybrid density functional calculations were done which suggests 
that the surfaces of small gold clusters are active cites for catalytic interactions. 
Among various systems they studied A1L13O2 clusters. In [22] electronic structure 
of 55 atom gold nanoparticles and catalytic activity from the modified electronic 
structure for small nanoclusters (~ 1.4 nm). In [23] density functional calculations 
of clusters were shown not to dissociate O2 

To study these effects we calculated the binding energy of the AU55O2 composite 
system using NWChem. The AU55O2 molecule is shown in Figure 7. The results 
for the binding potential are shown in Figure 8 and 9. These figures plot the energy 
as as a function of the distance of the center of mass of the O2 along the x axis 
using the coordinates of from Table 4. Figure 8 shows a very flat potential 

away from .8 nm so we show an expanded version of the potential in figure 9 from 
.805 to .845 nm. This shows a minimum at .815 nm and a weak binding energy of 
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Figure 7: A 55 atom gold icosahedron nanoparticle interacting with Oi- 




0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 
r [nm] 

Figure 8: Binding energy as a function of the distance to the center of mass of the 
O2 molecule to the center of mass of the AU55 nanoparticle along the x axis between 
.5 and .85 nm. 
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0.805 0.81 0.815 0.82 0.825 0.83 0.835 0.84 0.845 
Bond Length 

Figure 9: Binding energy as a function of the distance to the center of mass of the 
O2 molecule to the center of mass of the AU55 nanoparticle along the x axis between 
.805 and .845 nm. 
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-0.332135 eV. 

More interesting from the point of catalysis is to fix the center of mass distance 
of the O2, vary the distance between the two oxygens and study the energy as a 
function of this bond length. This is shown in Figure 10 together with the the 
energy as function of bond length for free standing O2 not in the presence of the 
gold nanoparticle. Here we fixed the distance from the origin to the center of mass 
of the O2 to be .815 nm. For all the data points that we studied in Figure 10 the 
presence of the Au^ nanoparticle lowered the energy and made the potential flatter 
as a function of separation of the two oxygens then would have been the case without 
the nanoparticle. This can be seen to give support to the idea that gold nanoparticle 
can be used to assist the breaking of O2, but from the figure the difference is small 
, less that .1 eV, so that more detailed studies are necessary to be definitive. One 
can also study the effect of different orientation of the O2 with respect to the gold 
nanoparticle as this can effect the energy of the system as well as the binding of the 
O2, rather than bring in the O2 along the x axis as we have done. 

5 Conclusion 

In this paper we have begun a study of the energy landscape for gold nanoclus- 
ters using density functional theory run on large parallel computers. By restricting 
ourselves to cuboctohedral and icosahedral geometries as well as intermediate ge- 
ometries we have explored a small subspace of low energy gold nanostructures. We 
found that the cuboctahedral shape was lower energy than the icosahedral shape for 
13 atoms and the opposite was the case for 55 and 147 atoms. We also estimated 
the energy difference. More significantly we were able to calculate the potential 
energy as a function of bond length and shape by introducing a shape parameter 
9. Finally we were able to calculate the effect of the presence of nanogold on the 
binding energy of O2 showing that it leads to weaker binding which may have ap- 
plications for the use of nanogold for catalysis. In the future it will be interesting 
to expand the configuration space of geometries that can be efficiently probed for 
gold nanoclusters. Because the chemistry and material properties of these clusters 
depend strongly on geometry this is likely to have practical applications. Rapidly 
expanding computational resources dedicated to nanoscience indicates that this will 
be possible in the near future. 
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